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O '. 1 Introduction 

(N 

■ This paper enters a work programme concerning the development of two-scale numerical 

Q ! methods to solve PDEs with oscillatory singular perturbations linked to physical phenom- 

ena. Recently, these methods have been tested on simple problems. For instance, in Ailliot, 
Frenod and Monbet[T], such a method is used to manage the tide oscillation for long term 
drift forecast of objects in the ocean; in Frenod, Salvarani and Sonnendriicker[5], such a 
method is used to simulate a charged particle beam in a periodic focusing channel. 



The ultimate goal of this work programme is to propose efficient numerical methods to 
simulate plasmas submitted to strong magnetic field. Of course, simulations of magnetic 
confinement fusion are part of this ultimate goal. 

Before going further towards this ultimate goal, several questions concerning the be- 
haviour of the concerned methods in front of non-linearities and non-smooth solutions need 
^ I to be investigated. The field of weakly compressible ID isentropic Euler equations offers a 

ly^ ■ relatively confortable framework in order to tackle these questions. 

G^ , Indeed, the mathematical context established in Grenier[9], Klainerman and MajdafTOl 

[TT] . Fortenbach, Frenod, Munz and Sonnendriicker [3], Frenod, Raviart and Sonnendriickerfl], 
Frenod and Sonnendriicker [H [7], Majda[14|, Metivier and Schochet[l6], Munz[l7], Schochet[2Q 
^ . EU [22j [22j , offers a nice marked out way. This is the first motivation of the present paper. 

I A second motivation originates from the fact that robust numerical methods such as 

^ ' finite volumes set out in Godunov[8j, LeVeque[l2] and Roe[l9] have been developped in 

order to solve compressible or incompressible Euler equations. Nevertheless, when they are 
^ , used to simulate the flow in an experiment with conditions inducing a small Mach number, 

?H I the CPU time cost of these methods is too high for operational applications because of the 

very small time step required in order to capture high speed travelling waves that appear 
in this case. 

The precise aim of this paper is to develop a two-scale numerical method to solve 
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where {x,t) £ x [0,T). In this model, e is the Mach number, u"^ = u'^{x,t) is the 
dimensionless speed of the fluid, = p^{x,t) is linked to the dimensionless density = 
p'^{x,t) by the relation 

and 

q'{p') = ^^{p'f + eqlip^), (1.3) 

where 7 is the adiabatic coefficient of the fluid and ql{p^) is regular. This form of the 
isentropic Euler equations is set out in Grenier[9] (p. 494) and Metivier and Schochet[16|. 
Equation (jl.ip is obtained from the classical form of the ID isentropic Euler equations 

dtu' + d,u' + ^d, {p{p')) = , 

dtp'' + da,{p'u') = , (1.4) 
ti^(x,0) = uo{x) , 
pHx,0) = l + epo{x) , 

where usually, pressure function is given by p{p'') = Using the relation p.2p between 

p*^ and p*^, we obtain 

Making and expansion of (1 + e p'^)^~'^, we obtain 



^- d^pip^)) = - (9,/) (1 + (7 - 2) ep^ + ql{p^)) , (1.6) 



g2 pe g 



where ql{p^) is regular, and, by introducing this result in the first equation of (|1.4|) . we 
deduce the first equation of (11.11) . In order to obtain the second equation of (II. ip . we 
introduce the relation (II. 2|) in the second equation of p.4p . 



In order to achieve this precise aim, inspired by two-scale convergence theory devel- 
opped in Nguetseng[T8] and Allaire|2j, we establish a homogenized model describing the 
mean behaviour of {u'',p''). This homogenized model neither contains nor generates high 
frequency oscillations but only their mean action. We also establish a way to reconstruct 
the oscillating solution {u'^,p'^) from obtained mean behaviour. The homogenized model 
and the reconstruction procedure involve an additional variable which is a copy of the time 
variable and which allows the transfert of informations between the mean behaviour and 
the high frequency oscillations. 



2 Results 

In this section, we claim the mathematical results leading to the construction of our two- 
scale numerical method. 

After recalling existence and properties of solutions of (|l.ll) . we give a first theorem 
giving the homogenized model and convergence properties. Then, we set out a numerical 
approximation for this homogenized model based on Roe's finite volume method which 
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is the two-scale numerical method that allows to solve Finally, we give a theorem 

giving some properties about the convergence of the considered numerical scheme. 



It is well known (see Grenier[9], Klainerman and Majda[l0], Majda[Tl], Metivier and 
SchochetpJl, or Schochet |20l [22l [23] ) that, with initial data satisfying 



there exists an existence time T > 0, independent of e, such that the system (II. ip admits 



aunique solution (u^p^) in (C(0, T; iJ^(T^)))^n (C^ (O, T; iJ^-i(Ti)))^ for any e > 0. Fur- 
thermore, the sequence (u^ p^)^ > o is bounded in (C {0,T; H%T^))fn {C^ (O, T; H'-^{T^))) ^ 



for the usual norm. 

Having this existence result at hand, we can claim the following theorem giving the 
homogenized model. 

Theorem 1. If we fix cq > and s > 2, there exists a constant K > such that, for any 
e G ]0, eo], we have 




(2.1) 




(2.2) 




where U and R are defined by 



U{x, T, t) = F{x -T,t) + B{x + r, t) + — 

zvr 

R{x, T, t) = F{x -T,t)- B{x + T,t) + -^ 



(2.3) 



with u = 




uq{x) dx, p 




po{x)dx, and F, S G C(0, T; i/*(T^)) satisfying 





(2.4) 



and 




(2.5) 
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A way to interpret this theorem is that if we can compute some good approximations 
Fh and of the solutions F and B of (|2.5p . we are able to reconstruct some functions Uh 
and Rh in the same way as we reconstruct U and R in (|2.3I1 . Then the obtained functions 
{x,t) ^ Uh{x,j,t) and {x,t) \—>- Rh{x,^,t) are good approximations of {x,t) f—>- u'^{x,t) 
and {x,t) 1-^ p'^{x,t) respectively. 

Based on this idea, we build our two-scale numerical method: firstly, we consider a 
uniform mesh on x [0, T] with space step h and time step k, and we use the notations 
Xi = ih and t„ = nk. Then we build F^ and B^ such that 



Fh{x,t) =F^ and Bh{x,t) = B^ V(x,t) G [xi_i/2, Xi+i/2[x , 
where and i?" are computed with Roe's finite volume approximation of (|2.5p : 



(2.6) 



7 + 1 

2/i i 



7 + 1 



+ 



4 

7 + 1 



+ Fn + 

(i^" + Fr_i) + 



47r 

2m + (7 - l)p 



2li + (7 - l)p 



4tt 



{F^ - FIU)\ , 



(2.7) 



Bl 



n+l 



n n2 /Rn , 2u - (7 - l)/9 



7 + 1 



+ 



(i^m + ^r) + 



^(i?r+i?r_,) + 



47r 

2n-(7-l)p 



47r 

2n-(7-l)p 



47r 



(2. 



Finally, we compute Uh and Rh by 



C/,,(x, r, t) = Fh{x -T,t) + Bh{x + T,t) + — , 

zvr 



Rhix, T, t) = Fh{x -T,t)- Bh{x + T,t) + 



27r' 



and we have the following convergence result. 



(2.9) 



Theorem 2. If s > 2, then the approximations Uh{-,T,-) and Rh{-,T,-) converge to 
U{-,T,-) and R{-,T,-) in L^{[0,T) x T^) for any r G Furthermore, if s > 3, the 
local truncation errors of the numerical scheme ^2.6\ }- ^27}\ }- /i2.8\} are first order accurate. 



Remark 3. The proof of the convergence result 12. S^) and of Theorem \E do not work if 
1^ 



s g]|,2[. However, in this case, the regularity is enough to get 



— U{-, -, •) ^0 weakly — * . 



R{; 



weakly — * . 



(2.10) 



m 



L°°(0,T;F^(Ti)) when e 0, with U and R solution of (E^j-flJ). 
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3 Construction of the two-scale numerical method 



3.1 Homogenized model 

We present here the construction of the model (12.5(1 . We recall that we start from the 
weakly compressible ID isentropic Euler equations (|l.ip - (|1.3p . Introducing the functions 
6*^ defined by 



nx,t) = ^(u%x+l,t)+p^{x+l,t)-^), 

b'{x,t) = Uu'{x + -,t) - p'{x + -,t) """^ 



(3.1) 

we can rewrite the system (11.11) : 

a,r{x, t) + 1 8x U {r(x, «) - i>'(i' + f ■ t) + ^ 

+ 8. {/■(.. 0) ^ + ^ (f + f . «)) ' + ^^^^^ /■(-.«) (3.2) 

8,b'(x, () + 5 8. US (/• (i - f , *) - ''"(^. *) + 1^) j 

+ 8il±i{6-(.,*))^ + _ + (3.3) 

equipped with 

/|:.o = /o = ^(^o + Po-^), (3.4) 

bU,=bo = l{uo-po-^). (3.5) 

Since the sequences (/^)e>o and {b'^)e>o are bounded in L°° (0,T; H^{T^)^ , there ex- 
ist two functions F and B in L°° (O, T; i?"(T^)) such that for some subsequences always 
denoted (/'^)e>o and (6'^)e>0; we have 

f^F weakly-* in L~(0,r;i7'(T^)) , (3.6) 

¥ weakly-* in (O, T; //"(T^)) , (3.7) 

for e — > 0. 

A first property of F and B is that their averages are equal to 0: to show this, we 
integrate (|3.2p - (|3.3p with respect to x to obtain 

dt(j^J'ix,t)dx] =dt( j^^b%x,t)dx] =0. 



Then we deduce that, for all t £ [0, T), we have 



/ fdx= fodx = 0, / b^dx= [ bodx = 0. (3.8) 

These results, combined with the fact that /"^ and b'^ weakly-* converge to F and B re- 
spectively, lead to the results (|2.4|) . 

Furthermore, usual compactness results (for example Lions[T3]) yield that the func- 
tional space 

U=[ge L-(0,T;i7^(Ti)) : dtg G {0,T; H'~\T^))} , (3.9) 

provided with the usual product norm, is compactly embedded in L°° (O, T; H'^^^{T^)) . As 
a consequence, since (/^)e>o and {b^)e>o are bounded in U, we have 

f'-^F strongly in L°°(0,r; //"-I (T^)) , (3.10) 



b' "-^ B strongly in L°°(0, T; //"-^(TI)) . (3.11) 

Since /"^ and ¥ are continuous in t, we can remark that F,B ^ C(0, T; //^^^(T^)) , so 
the convergence result (|3.10I) - (I3.11|) is also true in C(0, T; //''^^(T^)) . Having now this 
convergence result at hand, we can look for the constraint equations on F and B. For 
that, we multiply the equations (I3.2I )-( |X3]1 by a regular function ip with compact support 
on X [0, T) and we integrate with respect to x and t on x [0,T). We obtain 

/ f^{x,t) dt(p{x,t) dx dt — / fo{x) ip{x,0) dx 

JTi JTi 

I ^^{r{x,t))^dMx,t)dxdt 

^ ^7 ~ fi^^ t) dMx, t) dx dt 
47r 

+ dMx,t)dxdt (3.12) 





T 



Jt 

T 



4 V ^ e 



/" f '^—r^hUx + —,t)d^Lp{x,t)dxdt 

Jo Jti 47r e 

/ / ^^r(x,t)6^(x + — ,t)5^(/7(x,t)da;dt 

JO JTi ^ € 



dxdt = . 
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and 



¥ {x , t) dt^fix , t) dx dt — / bQ{x) ip{x,0) dx 
/ — - — [b'^{x,t)) dxip{x,t) dx dt 

2u - (7 - l)p t) dx dt 







47r 



JT 
T 
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(/'(a;-— dxip{x,t)dxdt (3.13) 



+ / / , r{x--,t)dM^,t)dxdt 
Jo Jji 47r e 
T 







^ f^ix — —,t) dxf{x, t) dx dt 

1 2 e 



J^^<ll{f'i''-J^i)-b"i^^t) + ^)d'M^,t)dxdt = 0. 

Because of (IS.lOj) and (|3.11l) . passing to the limit in the four first terms of (I3.12p and 
(|3.13p is straightforward. Because of the factor | in front of the last term of (|3.12|) and 
(I3.13p . we deduce that these terms converge to 0. 



of variables 



f f 2t 

In order to find the limit of / / b'^[x-\ , t) t) dt we firstly make a change 

3les 

/ / b^{x + —,t) dx(p{x,t)dxdt = [ [ ¥{x,t)dx 

Jo JTi £ Jo Jf^ 



(fix ,t] dxdt , 



then, we write 
rT 



f f 2t 

/ / b^[x-\ ,t) dxip{x,t) dx dt 

Jo Jti e 

fT r 2t f 2t 
= I I [b^[x,t) — B{x,t)) dx(p{x ,t)dxdt+ / / B{x,t)dx(p{x ,t^dxdt. 

Then, using the Cauchy-Schwarz inequality, we show that 

f 2t 

lim / / {b^{x,t) - B{x,t)) dx^{x ,t)dxdt = 0. (3.14) 



If we define ip'^{x,t) = (p{x,j,t) where (p{x,T,t) = (p{x — 2r,t), the sequence {cp'^)e>o is 
bounded in L~(0, T; H'{T^)). Then, there exists a function $ in (O, T; L'^{T^; H'{T^))) 
such that 

lim / / 'ilj[x,-,t) ip^{x,t)dxdt = / / / ip{x,T,t) ^{x,T,t) dr dx dt (3.15) 
<^^oJo Jt^ f Jo Jti Jti 

for any regular function ^p defined on x x [0, T). As described in Allaire[2], $ is 
called the two-scale limit of f^. Furthermore, following Marusic-Paloka and Piatnitski[15| 
and Allaire[2], since if is regular and 27r-periodic in r, it is an easy game to prove that 
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^ = (pin L°°{0,T;L^{T^;H'{T^))). As a consequence, we have ^{x,T,t) = ip{x - 2r,i) 
and 

(p' ^ [ ip{x - 2T,t)dT weakly-* in L°°{0,T;H'{T^)) . (3.16) 
Using these results, we obtain 

B{x, t)dxp{x ,t) dx dt 

= -/ / / {dj:B{x,t))ip{x-2T,t)dTdxdt (3.17) 

Jo Jt^ JTi 

= -\j^ j B{x,t){^j dr{(pix -2T,t)) drj dxdt = 0, 

and finally 

n2t 
b'(x + —,t)da,ip{x,t)dxdt = 0. (3.18) 
_'i e ' 

In the same way, we obtain 

n2t 
f(x ,t)dxip{x,t)dxdt = 0, (3.19) 

lim / / (b^{x+—,t)) dxp{x,t)dxdt = 0, (3.20) 
rT 



e->0 




2t ,\2 



-'Ti 



lim / / lf^{x ,t)) dxip{x,t)dxdt = 0. (3.21) 



n2t 
f^{x, t)b^(^x -\ , t) dx<p{x, t) dx dt, we notice that 

n2t 
f{x,t) b^{x + —,t) dxip{x,t)dxdt 

= / / {f{x,t)-F{x,t))B{x+—,t)dx(p{x,t)dxdt 

Jo JTi C 

n/ 2t 2t \ 

F{x,t) [b^{x + —,t)- B{x + —,t)j dx(p{x,t)dxdt 



j-T 

+ 



n2t 
F{x, t)B{x + —,t) dxpix, t) dx dt . 



With Cauchy-Schwarz and Holder's inequalities, we easily show that 

n{rix,t) - F{x,t)) B{x + —,t) dxp{x,t)dxdt ^0, (3.22) 

n/ 2t 2t \ 

F{x,t) {b^{x + —,t) - B{x+ —,t)) dxip{x,t)dxdt ^0, (3.23) 
■1 V e e / 

J {f{x,t)-F{x,t))(b''{x+'^,t)-B{x + ^,t)^dx(p{x,t)dxdt^O, (3.24) 



when e — > 0. If we define the functions B^{x,t) = B[x + Y'*)) remark that the 
sequence (-B'^)e>o is bounded in L°° {0,T; H'^^^ (T^)) and then admits a two-scale limit 5B 
in L'^(0,r;L°°(Ti;F^-i(Ti))) satisfying 

^{x,T,t) = B{x + 2T,t) V(x,r,t) G X X [0,T) , (3.25) 

^ [ B{x + 2T,t)dT = weakly-* in L°°(0,r;i7"(T^)) . (3.26) 
This means that, for any function ^p £ C{0,T; H^^^(T^)^ , we have 

lim / / B'{x,t)'ip{x,t)dxdt = 0. (3.27) 
Setting then ■0 = F dx^ in (I3.27p . we obtain 

hm / / F{x,t)B{x + —,t)dx'^{x,t)dxdt = 0. (3.28) 



Hence 



/" /■ 2t 
lim/ / /'(x,t)6'(x + — ,t) a^.(/9(x,t)dxdt = 0. (3.29) 



The same method also gives 



lim /" / b^{x,t) r(x -—,t) dx(p{x,t)dxdt = 0. (3.30) 
With these results, we can pass to the limit in (I3.12|) - (I3.13I) and get 

F{x,t) dt(p{x,t) dx dt — / fo{x)ip{x,0)dx 
I' I, {^{F{x,t)f + ^^^^^^^^^^F{x,t))dM^,t)dxdt = 0., 



T 



JT'^ 
T 



and 

r-T 



[ / B{x,t) dt(p{x,t) dx dt — / bo{x) ip{x , 0) dx 

iTi JTi 

r I {B{x,t)? + ^^^-^^^^^^ B{x,t)) dMx,t)dxdt = . 

Jo Jti ^ 4 ' 47r / 



(3.31) 



(3.32) 



Remembering the definition (I3.4I) - (I3.5I) of /o and bo, we recognize here the weak formulation 
of (1231) . 

Under the hypothesis of Theorem [1] about the initial data in (12.51) . we deduce that the 
solution {F,B) belongs to {C{0,T; H''{T^))) and is unique in this space (we can apply 
the Theorem 3.6.1 of Serre[24j for example). Finally, from uniqueness, we have all the 
convergence results above for the whole sequences {f^)e>o, ib^)e>o and not only for some 
subsequences. 
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3.2 Convergence of the homogenized model: proof of Theorem [T] 



In this section, we will prove inequalities (12. 2p . Before going further, we remark that con- 
vergence (I3.6I) - (I3.7I1 is already proved and that, thanks to (12.31) and (|3.1I) . these inequalities 
are equivalent to the theorem below. 

Theorem 4. For eo > and s > 2 fixed, there exists a constant K > such that 

\\n;t)-Fi;t)\\l,^^,^dt^ <Ke, (3.33) 

p{;t)-B{;t)\\l,^^,^dt^ <Ke, (3.34) 

for any e g]0, eo]. 

Proof of theorem^- to prove this theorem, in the same spirit of Frenod, Raviart and 
Sonnendrucker|4j, we introduce the functions 7*^ and 6^ defined by 



Y{x,t) = -^{r{x,t) - F{x,t)) - W{x,*-,t) , 
6'{x,t) = ^{B{x,t) -b'{x,t)) -V{x,*-,t) , 



(3.35) 



with W and V defined on x x [0, T) by 

W{x,r,t) = -^^^ J d^[{B{x + 2e,t)f -2F{x,t)B{x + 2e,t) 



-B{x + 2e,t)]d0. 

TT 



V{x,T,t) = ^^j^ d^[{F{x -2e,t)f -2F{x -2e,t)B{x,t) 



(3.36) 



+ -F{x-2e,t)]de . 

TT 



Concerning sequences (7'^)eg]o,eo] and (<^'^)ee]o,eo]' have the property below. 

Lemma 5. Under the hypothesis of theorem^ there exists a contant M > 0, independent 
of e £ ]0, eo] and t £ [0, T), such that 



^*(l|7'(-,i)||i2(Tri) + ||'^'(-,*)|li2(Tri)) 



(3.37) 

<M(^ + \/2](||7^(.,i)||^,(^,) + ||<^^(.,t)||^,(^,)) +^/2M, 

for any e G ]0, eo] and t G [0, T) . 

Having this inequality at hand, we can apply Gronwall's lemma to find a constant 
L > 0, independent of t and e, such that 

\Hi-M%iT^) + ¥i-M%(T^)<L^ V(*'e) G [0,r)x]0,eo]. (3.38) 
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Finally, we define the constant K by 



K 



te [o,T) 

r eTl 



i G [0,T) 



(3.39) 



and we obtain the inequalities (|3.33p - (|3.34p . giving the theorem. 

Now, to finish the proof of Theorem |4j we have to prove Lemma [5l Firstly, we notice 
that dt'y'^ and dtS'^ read 



2t \ 

a{r{x,t) + F{x,t)) + CB{x + -,t) -p+)d^Y{x,t) 

2t 

€ 



+ [ - ad4r{x,t) + F{x,t)) + Cd^B{x + j,t)j j'{x,t) 



+ UB{x + -,t)+b^{x + -,t)- 2 fix, t) - dj^ {x + -,t) (3-40) 
2\ e e IT/ e 

t / 2f 2f \ 2f 

+ I 9, + -, t) +b'{x + -,t) - 2r{x, t)) 5'{x + -,t) 

+ T'{x,t), 



/ 2t \ 

dtS'ix, = ( - a{B{x, t) + t)) + CF{x - —,t) - /?_ j t) 

+ (^-ad.4B{x,t) + b'{x,t))+CdxF{x-^,t)')6'{x,t) 

+ (x - - + F(x - - 26^(x,t) + I) d.Y{x --,t) (3.41) 

+ ^ 5. (r (x - ^, t) + F(x - |, t) - 26^(x, t)) - f , 
+ A^(x,t), 

with a, C, r'^(x, t) and A'^{x,t) defined by 

^.-^^i^. C^^. (3.42, 



r^(x, t)=- ad, (W{x, -,t) (r (x, t) + F(x, t))) - p+d,w{x, -,t) 

- Cda^[nx,t)V{x + ^,*-,t) - B{x + ^,t)W{x,*-,t)) 

+ 1 a. (^y(x + 1, ^, t) (s(x + 1, t) + b%x + 1, t))^ 

- g (x + |, ^, t) - i 9 JgJ (/ (x + ^, 0) ) - (x, l,t) 
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and 



A'{x, t)=- ad:, {V{x, -,t) {B{x, t) + ¥{x, t)) - /3-d^V{x, -,t) 

- Cd^(w{x -j,^-,t)b^{x,t) ~ F{x -j,t)V{x,^-,t)) 

44(H-(.-f4.o(r(-f.0+n-f..))) <^"^' 

+ g 8,w{x - ^, + 1 8. L (p.(j, _ i , ()) j _ (a, I/) £, t) . 

Since the sequences (/"^)ee]o,eo] ^^^1 (6'^)£g]o,eQ] are bounded in C[0,T; H^(T^)^ with 
s > 2, F and S are in C(0, T; iJ*(T^)) , and is regular, we can find some constants 
Cj > (i = 1, . . . , 8), independent of e and t, such that 

2 / ( - a{r{x,t) + F{x,t)) - (5+ + CB{x + -,t)) Y{x,t)d^Yix,t)dx 
7ti ^ e y (3 45) 

< C'i||7^(->*)|L2(Ti) > 

2 /" ( -a(b'{x,t)+B{x,t)) - p^ + CF(x-—,t))s'{x,t)dJ'{x,t)dx 

<^2||<^^(-,i)||i.(T.), 

2 / (-aa4r(x,t) + F(x,t)) +C5a.5(x + -,i)) (7^(x,t))'dx 
Jti V ^ e ^ (3.47) 

<^3||7^(-,0||V). 



<C4\\n,t)\\%^^,^, 



(3.48) 



2 / + + -2f (x,t)) 5^(x + -,t)7^(x,t)dx 

7ti 2 V ' e e ^ e (3.49) 

<^5||7^(-'*)IL^(Ti)IK'(-'*)|L^(Ti)' 

f r / 2f 2/ N 2/ 

2 / iajF(x--,i)+r(x--,t)-26^(x,0)7l^--,i)<^'(^,0rf^ 

2 V ^ e ^ ^ e ^ J ^ € ' (3.50) 

<C6||7^(-,0IL.(P)II'^^(->*)IL^(T1)' 

2 j^r{x,t)^'{x,t)dx < C7||7'(-,0|L2(Ti) , (3.51) 

2^^ A'{x,t)S'{x,t)dx < Cs\\S'{;t)\\^,^^,^ . (3.52) 

On the other hand, making in the first integral the change of variable x ^ x + t and in 
the second one x i— > x — r, rearranging the terms and using the definitions of 7^ and S^, 
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we have 



C [ (b'{x + 2r, t) + B{x + 2r, t) - 2f{x, t) - -) dj'{x + 2r, t) 7'(x, 

J^i V vr/ 

+ C / fr(a;-2r,t) + F(x-2r,t)-26^(x,t) + ^)9a;7'(x-2r,t)(5'(x,t)dx 
= C / (6'(x + r,i) -S(x + T,i))7^(x-r,t)aa;5'(x + r,t)<ix 
+ 2C / - r, - /'(x - r, t))7'(x - r, t) dj^{x + r, dx 

-2C / fF(x-r,i)--B(x + r,i) + — ')7'(a;-r,t)5a;5'(x + r,i)dx 
+ C / (/'(x - r, - F{x - T, t))S^{x + T, t) d^cl^ix - T, t) dx 
+ 2C / {B{x + T,t)-b'{x + T,t))6'{x + T,t)d^j'{x- T, t) dx 
+ 2C / (f(x -T,t)- B{x + T,t) + ^) 6'{x + T, t) d^Y{x - T, t) dx 
= C / e(5x7'(a;-T,t)-aj;(5'(a; + r,t))7'(x-r,t)(5'(x + r,t)dx 
+ C / e5a;7'(x - r, t) {W{x - r, r, t) + 2F(x + r, r, i)) + r, t) dx 
-C / eaa;(5'(x + r,t)(y(a; + r,r,t) + 2PF(a;-r,r,t))7'(x-r,t)dx 
-2C / ea^(5"(x + T,t)(7"(x-T,t))^(ia; 
+ 2C / €d^-f'{x - T,t){5'{x + T,t))'^ dx 

+ 2C / r, t) {S^{x + r, t)0^7'(x -T,t)- 7'(x - r, t)dj''{x + r, t)) dx . 

The sequences (€7*^)^ gjo^eg] and (65*^)^ gjOjCo] ^-r^ bounded in C(0, T; if*(T^)) , so we can find 
some constants Cj > (i = 9, . . . , 13), independent of e and t, such that 

C / ({d:cl%X-T,t)-dJ%X + T,t))-i'^{x-T,t)5%X + T,t)dx 

< C9||7^(-, ^2(Ti) l|7'(-, *)|| i2(Ti) , 

C / ed:c7^{x-T,t){W{x-T,T,t) + 2V{x + T,T,t))6^{x + T,t)dx 

<^^io||<^^(-,i)|L.(^,), 

C / edJ^{x + T,t){V{x + T,T,t) + 2W{x-T,T,t))-f^{x-T,t)dx 

<Cn||7^(-,i)||^,(^,^, 

2C e5^(5^ (x + T, t) (7^ (x - r, t)) ' dx < C12 1 1 7' (• , 1 1 5,2 (Ti) , (3-56) 



(3.53) 



(3.54) 



(3.55) 
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2C l^^ed,Y{x-T,t){5'{x + T,t))''dx<Cis\\5'{;t)\\l,^^,y (3.57) 
For the last integral, we simply compute: 

2C / R{x, T, t) {5^{x + r, t)d^Yix - T,t) - Yix - r, t)dx6'{x + r, t)) dx 

= 2C / dr {R{x, T, t)) 6^ {x + T, t) Y {x - T, t) dx 

-2(5^1 / R{x,T,t)6^{x + T,t)Yix-T,t)dx] 

\Jt^ J (3.58) 

= 2C / dr{R{x,T,t))6'{x + T,t) 7*^(3; — r, t) dx 

= -2C / S^{x + r, t) 7'(x - r, t) ^^^^/(a;, r, t) dx 

< C'u ||7^(-)*) II i2(Ti)ll'^^(''^)IL2(']ri) 

for a constant C14 > independent of e and i. 

Having (l3^ - ([338]) at hand, we multiply ([3^ by -f'{x,t) and (l3^ by (5'(x,t) to 
deduce 



5, 



'*(lb^(-'*)|lL^(Ti) + ll^^(-'*)|lL(Ti)) < ^(l|7^(-'*)|li.(Ti) + II^^(-'*)|Il^(T1) 

+ l|7^(-'0|| L2(Ti) + ll'^^("'*)|L2(Ti) (3-59) 
"l~ II7 ('' ^) llia^xi) 11^ ^)|Il2('i['1)) ' 

with M = max( Ci + C3 + C12, C2 + C4 + Ci3,C 5 + Q + Cg + Cu, C7 + Cn, Cs + Cio). 
Inequality (|3.37p is directly obtained from (|3.59|) . ending the proof of lemma [5] and then 
the proof of theorem HI □ 



3.3 Properties of the finite volume scheme: proof of the Theorem 2 

Thanks to ([231) and (I2J|), convergence of {Uh{-,T, ■), Rhi-,T,t)) to ([/(•, r, •), -R(-, t, t)) for 
any r G is equivalent to the convergence of {Ffi, B^) to (F, B). 
Furthermore, the equations satisfied by F and B are of the form 

dtq + d^{f{q)) = q{x,^)=qQ{x), [ qix,t)dx = 0, Vt, (3.60) 

with f{q) = aq^ + Pq, a and /3 = P± defined in (|3.42p . In the same way, the numerical 
method (1221)- ([221) "dSISl) is of the form 

Qh{x,t) = Qi y[Xi^i/2,Xi^i/2[x[tn,tn+l[, 

Qr' = Q1- l{HQ7+i,Q7) - HQf,Q7-i)) , (3.61) 
HQi, Q^-l) = lifiQi) + f{Q^-l)) - \HQi + g.-i) + p\{Qi - Qi-i) , 
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27r 

with a space step h = — and a time step k. 

We notice that the CFL condition at the n-th time step is given by 

^ max |a(gr + Qti) + P\ = < I ■ (3.62) 

With these notations, proving Theorem 2 is equivalent to proving the theorem below. 

Theorem 6. If s > 2, the approximation Qh converges to the solution q of 113. 60\) in 
-L^([0,T) X T^) norm. Furthermore, if s > 3, the local truncation error of the numerical 
scheme 113. 61\) is first order accurate. 



Proof of theorem 5: the total variation of a function q £ L^([0, T) x T^) is defined by 
v^o V Jo Jti l' 

For an approximation Qh computed with the numerical method (I3.6ip . we have 



TVt (q) = lim. sup - f f \q{x + r],t) — q{x,t)\ + \q{x,t + rj) — q{x,t)\ dxdt. (3.63) 
v^o V Jo Jt^ ^ -I 



T/k 

TVriQh) = Yl [kTV{Q'') + - Q"||i] , (3.64) 

n = 

with TV{Q) and ||Q||i defined by 

Nx Nx 

TViQ) = Y,\Qi-Qi~i\ and ||Q||i = /i ^ |Qi| • (3.65) 

1=0 1=0 

Introducing the sets £. and /C defined by 

£ = |g G L^([0,T) X T^) : j gdx = OVt}, IC = {q e £ : TVr{q) < R} , (3.66) 
where i? is a constant depending on qo, it is well known that /C is a compact subset of C. 

Since the numerical flux involved in (|3.61|1 is continuous and satisfies 

HQ,Q) = f{Q) Vg, (3.67) 

the considered scheme is consistent with the conservation law (I3.60p . Hence, applying 
LeVeque[l2], proving that Qh converges to the solution q of (|3.60ll in L^([0,T) x T^) norm 
for /i — > reduces to prove that the scheme is TV-stable. In other words, we need to prove 
that for h G [0, ho] , the approximation Qh lies in some fixed set /C where R only depends 
on the initial data qo, the final time T and the function /. 

The TV-stability is the consequence of two lemmas. 
Lemma 7. The numerical method i3. 61\) is TVD (Total Variation Diminishing), i.e. 

ry(g"+i) < ry(Q") Vn. (3.68) 
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Proof of the lemma\^ we develop Q^^^ - Q^: 

QUI - QU' = - Q7 - |^(/(^?r+2) - fm - /(Qr+i) + /(Qr^i) 

+2|a(Q[Vi+Qn+/?|(Qm-Qn 
-|a(Qr+Q?_i)+/?|(Q?-Qr-i)) 

= (i - sg(«(Qr+2 + Qr+i) + («(Qr+2 + Qr+i) + p) (Qr+2 - ^r+i) 
+ 1^ (i + sg(a(gr + Qr-i) + /S)) (a(Qr + Qr-i) + /?) (q? - qu) , 

where sg stands for the usual sign function. Using the CFL condition (|3.62|) . we write: 

\QUl - QU'\ < 4 (i - sg(«(Qr+2 + Qi+i) + /?)) |"(Q?+2 + + /3| IQIV2 - Qr+ii 
+ (1 - ^ |a(Q[Vi + + iQr+i - Qri 

We deduce then 

TViQ-+') = ^ \QU! - Qr^l < E l^^i - = ^^(^") (3.69) 
giving the lemma. □ 

As a consequence of this lemma, we deduce that there exists a constant Mi > which 
only depends on the initial data go and such that 

TV{Q'')<Mi Vn. (3.70) 
Lemma 8. For any n, we have \\Q^'^'^ — Q"||i < 2Mih. 

Proof of lemmalM' we develop ~ Qi- 

- f{Q7) + f{Q7-i) + HQ7 + QU) + /?! {Ql - QU)) 

= (1 - sg(a(Qr+i + Q7) + /?)) + QD + /3) (Qr+i - qd 
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Then, using the CFL condition, we obtain 



|Q«+l-Q"||^ = /i^|Q^+l-Q^| 



i = 



< 



i = 

i = 

<2hY,\Q7-Q7-i\ 

i = 

and we complete the proof by applying lemma [71 □ 



-|a(Qr+i + Qr)+/3|IQm-Q? 



-\a{Q2 + Qti) + P\\Q?-Q7-i\ 



End of the proof of theorem 5: combining all these results, we finally have the inequality 

TFr(Qh) < (l + ^)rMi. (3.71) 

Remarking that ^ is bounded by a constant M2 which only depends on the initial data go, 
we finally obtain that there exists a constant C > which only depends on go and T, and 
such that 

TVTiQh) < C (3.72) 

for any space step h. We conclude that the numerical method (|3.6ip is TV-stable, yielding 
the convergence of Qh to q. 



Now, we have to prove that, assuming that s > 3, the local truncation error is first 
order accurate. With the notation above, the expression of the local truncation error reads 

n _ n+1 1 

+ |a(gr+ei)+/3|(9r-ei)), 

where qf = q{xi, tn)- According to the signs of a{qf + qf_i) + P and a(gJYi + + 
distinguish four different cases. 



First case: we assume that a(g"+i + qf) + /? > and a{q'^ + qf_i) + /3 > 0. Then, 
the local truncation error becomes 

.(..,,„)^gL^- /W)-/W-.) , ,3.,,, 

Since s > 3, g G C^([0,T) x T-*^), so we can write 

e(x„t„) = -dtq{xi,tn)-dJf(qixi,tn)))+Oik)+Oih)=Oih + k). (3.75) 
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Second case: we assume that + qf) + /? < and a{qf + qf_i) + (3 <0. Then, 

the local truncation error becomes 

Since s > 3, g G C^([0, T) x T^), so we can write 

e(x„ t„) = -dtq{xi,tn) - (f{qixi, t„))) + Oik) + = + k) . (3.77) 



• Third case: we assume that a((?"+i + qf) + P > and a(g" + (/j'L]^) + /3 < 0. Then, 
the function A defined by 

X{x) = a{q{x,tn)+q{x-h,tn)) +P (3.78) 

admits a zero denoted in [xi, Xj+i]. Since s > 3, A is of class and verifies 

X{xi) = X{x^) + 0{xi - x^) = 0{h) , ^^^^^ 
A(xi+i) = A(x*) + 0{xi+i - X*) = 0{h) . 

In the same way, we have 

/'(qD = 2aq? + P = A(x,+i) + Oih) = 0{h) . (3.80) 

The local truncation error is of the form 

n _ n+1 

eixi, tn) = = -dtqix^,tn) + Oik) = f'{qix^,tn)) d^qix^,tn) + Oik) . 

(3.81) 

Then, using (I3.80p . we obtain 

eixi,tn) = Oih + k). (3.82) 

• Fourth case: we assume that a(g"+i + qf) + P < and aiqf + q'^_i) + jS > 0. We 
proceed as we did in the precedent case, i.e. we deduce that the function A defined 
by (I3.78P admits a zero x* G [xj,Xj+i], and then the result (I3.80p is true. Hence, 

e{Xi,tn) - ^ ^ 

= -dtqixi,tn) -2d^(^f{qixi,tn))^ +Oih + k) (3 33) 

= -/'fe")9.(?(xi,t„) + 0(/i + A;) 
= Oih + k). 

We conclude that, in any case, the local truncation error is Oih + k), i.e. first order 
accurate. □ 



4 Numerical results 

The first goal of this section is to numerically show that u'^ix,t) — Uhix, = ©(e) and 
p'^ix, t) — Rhix, ^,t) = 0(e). Secondly, we will briefly analyze the gain in term of CPU time 
our two-scale numerical method brings when compared with a classical method consisting 
in solving directly p.ip . Lastly, by using our method, we will explore the simulation of 
experiments which were not accessible before because of too small Mach number. 
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4.1 Convergence in e 



We consider here the initial data uq and po defined by 



, , , cos(x) , / X -, sin(x) 
uo{x) = 1 + and po{x) = 1 + 



(4.1) 



and we take a uniform mesh with + 1 points denoted xq, . . . ,xn^. 

As illustrated in figures 1 and 2, the functions Uh{x, ^,t) and Rh{x, ^,t) are very close 
to and respectively. The numerical experiment showed in these figures is made with 
e = 0.05, 7=1, and = 1023, T = 2.5. 




2.4?7468 and epsilon = 0.05 




0.90 1.79 2.69 



5.3K 6.2! 



0.00 0.90 1.79 2.69 3.59 



5.3S b.2l 



Figure 1: Comparison of u'^{-,t) and Uh{-, i,t) at times t = and t = 2.5 



I = 2.4y74M and epsilon = 0.05 




Z69 3.59 




2.69 3.59 



Figure 2: Comparison of p^{--,t) and Rh{--, \,t) at times t = and t = 2.5 



In order to quantify this good accuracy of Uh{x, j,t) with u'^ and Rh{x, j,t) with p*^, we 
compute the errors u'^{x,t) — Uh{x, and p'^{x,t) — Rh{x, ^,t) in LP([0, T) x T-*^) norm 
(p = 1,2,00), for several values of e ranging in {0.01,0.03,0.05,0.07,0.1}, and we obtain 
the figures and the array below: 
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0.40 



0.36- 
0.32- 
0.28- 
0.24- 
0.20- 
0.16- 
0.12- 
0.08- 
0.04- 
— 



0.08 



0.10 



O error in L'^l norm 
+ error in L'^2 norm 
X error in Max norm 



K»cpsilon with K = 2.8893757 
K»epsilon with K = 1.0358362 
K'epsilon with K = 1.3792216 



Figure 3: Error u'^{x,t) — Uh{x, ^,t) in L^, and L°° norms. 



0.40 
0.36 





















































































0.02 

O error in L'^1 norm 
+ error in L'^2 norm 
X error in Max norm 



— K»epsilon with K = 3.5842834 

— K»epsilon with K = 1.1780144 

— K'epsilon with K = 1 .690496 



Figure 4: Error p^{x,t) — Rh{x, j,t) in L^, and L°° norms. 
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Value 
of e 


Error u'^{x, t) — Uh{x, |, t) 


Error p''{x,t) - Rh{x, 


norm 


L"^ norm 


norm 


norm 


L'^ norm 


norm 


0.1 


0.2880016 


0.1034517 


0.1201132 


0.3634635 


0.1192177 


0.1651374 


0.07 


0.2043288 


0.0735498 


0.1160378 


0.2516643 


0.0828920 


0.1203342 


0.05 


0.1452756 


0.0518309 


0.0719549 


0.1744198 


0.0575395 


0.0879887 


0.03 


0.0845621 


0.0296311 


0.0483544 


0.0987818 


0.0326775 


0.0499608 


0.01 


0.0260695 


0.0085251 


0.0195340 


0.0303524 


0.0092954 


0.0269734 



Table 1: Errors u''{x,t) - Uh{x,\,t) and p'^{x,t) - Rh{x,\,t). 

The results above show an error u'^{x,t) — Uh{x, ^,t) in norm decreasing when e — > 
as Kie with Ki w 2.8893757, as with K2 ^ 1.0358362 in L'^{[0,T) x T^) norm, 
and and as K^e with K^o ^ 1.3792216 in L°°([0,r) x T^) norm. Concerning the error 
p'^{x,t) — Rh{x, ^,t), we observe that it decreases when e — > as K(e with K[ « 3.5842834 
in L^{[0,T) X norm, as K!^e with K!^ ^ 1.1780144 in L'^{[0,T) x T^) norm, and and 
as K'^e with K'^ w 1.690496 in L°°([0,T) x T^) norm. 



4.2 CPU time cost 

One of the motivations of the present paper is the high CPU time cost of classical methods 
when used with small Mach number, as we explained it in the introduction. For example, if 
we apply Roe's finite volume scheme on the model p.ip . we must verify the CFL condition 



k 

— max 

h i ,n 



< 1 , (4.2) 



where h is the space step, k is the time step, and Pj are defined by 



'v'r+W+<-iA/r+^ 



,/T+7pf+^l + epU 

= \ 7e{pr-Pr-i) ^ ^^-1 ' 

' [ (l + epf)^-i else. 

Hence, in order to garantee the stability of the scheme, the smaller e is, the smaller 
k must be (and the higher the CPU time cost will be). One advantage of our two-scale 
numerical method is that its CFL condition ()3.62|) does not depend on e, so k does not 
have to diminish with e. 



Value 
of e 


Roe's method 
on ([11]) 


Two-scale numerical method on (|2.3I)-(I2.5I) 


Computation of {Fh,Bh) 


Computation of {Uh,Rh) 


0.1 


24 m 34 s 36" 


3 m 45 s 31" 


42 s 79" 


0.07 


32 m 55 s 44" 


3 m 45 s 22" 


42 s 93" 


0.05 


44 m 7 s 21" 


3 m 45 s 38" 


42 s 74" 


0.03 


1 h 10 m 15 s 


3 m 45 s 34" 


42 s 80" 


0.01 


3 h 20 m 53 s 78" 


3 m 45 s 12" 


42 s 91" 



Table 2: Comparison in terms of CPU time cost between Roe's method on (jl.ip and the 
two-scale numerical method on (|2.3p - (|2.5p with the initial data (|4.1|) . = 1023, 7 = 1. 
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The array above gives the comparison in terms of CPU time cost between Roe's method 
apphed on the non-homogenized model (jl.ip . and our two-scale numerical method. These 
results have been obtained with the computer characteristics below: 

• processor: Pentium® M 715, 

• memory: 512 Mo DDR-RAM, 

• operating system: SuSe® 9.1 Pro, 

• compiler: gcc 3.3.3-41. 

In particular, we distinguish the time cost for the computation of i<X and on one 
hand, and the time cost for the reconstruction of Uh{x, and Rh{x, ^,t) on the other 
hand. 



4.3 Numerical experiments with small Mach number 

In this paragraph, we present some numerical results obtained with our two-scale numerical 
method on simulations of experiments inducing a very small Mach number. These results 
are obtained with the initial data (|4.1|) . + 1 = 1024 points in and 7=1. 





OUOU 0.62S 1.355 1.8S.1 1511 3.139 3.766 4.394 5.022 5.649 62 0.000 62S 1.255 I.SS3 2.511 3.1.19 3.766 4.394 5.022 5 64') 6.277 



Figure 5: Uhix, (left) and Rh{x, f ,t) (right) with e = 10"^ at times t = 0,t = l. 
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Function x -> R|x,l 





062S 1.255 I.8S3 2.511 3.139 3.766 4.394 5.032 5.649 6.2 000 0.628 1.255 1.883 2.511 3.139 3.766 4.394 5.022 5.649 6.277 



Figure 6: Uh{x, ^,t) (left) and Rh{x, (right) with e = 10"^ at time t = 2. 



Fimciioii X -> U(x,r/BpsiloD,i) with t = and epsih 



Fimctjod X -> R(ift/fipsiloD,t) with 1 = and epsiloa = 





FuEiction X -> U(x,r/fipsiloD,i) with t = 1 and epsiloD 



Functioii X -> R(x,t/fipsilon,t) with 1=1 and epsiloii = 




2.511 3.139 3.766 4.394 5.022 5.649 6.277 



Function x -> U(i.t/qisiloo,t) with t = 2 and qisiloii = 



Function x -> R(x^t/epsilon,t) vji^ t = 2 and epsilon = 




2.511 3.139 3.766 4.394 5.022 5.649 6.2 




Figure 7: Uh{x, l,t) (left) and Rh{x, l,t) (right) with e = lO"^ at times t = 0,t = l,t = 2. 



23 



Punedon i -> U(K,Uepailoii,i) with t = and epsilon = 0.000001 




(!.()«) ()h2N 1.255 I.SS3 2311 3.139 3.766 4.394 5.022 5.649 6.2 0.000 0.628 1.255 1.883 2.511 3.139 3.766 4.394 5.022 5.649 6,277 



PuDcdod I -> U(s,t^ep9iloii,i) widi t = 1 and epsilon = 0.000001 Punedon i -s HJitt/fpsilon,!) widi t = 1 and epsdon = 0.000001 




0.000 0.628 1.255 1.S83 2311 3.139 3.766 4.394 5.022 5.649 6.2 0.000 0.628 1.255 1.883 2.511 3.139 3.766 4,394 5.022 5.649 6,277 
U(x,l/epsUon,l) mtt/epsdon,!) 



Punedon i -> U(x.t/epailon.t) widi t = 2 and epsilon = 0.000001 Punedon i -> R(i,t/r^silon,i) widi t = 2 and epsUon = 0.000001 




0.000 0.62K 1.255 1.SS3 2.511 3.139 3.766 4.394 5.022 5.649 6.2 0.000 0.628 1.255 1.883 2.511 3.139 3.766 4.394 5.022 5.649 6.2T7 
U(x,l/epsilon,l) R(ii,t/(psilon,t) 



Figure 8: Uh{x, f ,t) (left) and Rh{x, f ,t) (right) with e = lO"*^ at times t = 0,t = l,t = 2. 

We can notice that these results have been obtained after a few minutes of computation 
just as described in the table below: 



Value of e 


Computation of (Fh.Bf,) 


Computation of (Uh, Rh) 


10- i 


3 m 4-5 s 71" 


42 s 84" 




3 m 45 s 95" 


42 s 87" 


10"^ 


3 m 45 s 75" 


42 s 98" 



Table 3: CPU time costs. 
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Furthermore, if we continue the simulation to final time T ~ 3.2 for e ranging in 
{10~\ 10-2, 10-^ 10"^, 10-^ 10-6}, we remark that discontinuities appear at time t ^ 
3.11 independently of e and that, for e ranging in {0.01,0.03,0.05,0.07,0.1}, these discon- 
tinuities also appear in the solution {u^,p^) of p.ip at the same time. 

5 Conclusion 

In this paper, we have developped a two-scale numerical method for the weakly compress- 
ible ID Euler equations. By using two-scale convergence tools from Nguetseng|18| and 
Allairep], we have developped a new model independent of the Mach number e. Then, 
using a finite volume scheme on this model, we have obtained numerical approximations 
that converge to the solutions of the weakly compressible ID Euler equation when we refine 
the mesh on T-*^ x [0, T] and when e ^ 0. 

Furthermore, we have proved that this two-scale numerical method is first order accu- 
rate in e. One the other hand, numerical results have confirmed this accuracy in e and also 
have proved that, even with a very small Mach number, the two-scale numerical method 
allows us to simulate experiments with a reasonable CPU time cost. 

Motived by the behaviour of the two-scale method, we are now ready to develop such 
a method on a model which describes a plasma submitted to a strong magnetic field such 
as the Vlasov-Poisson model established in Frenod and Sonnendriicker [7j. Furthermore, 
we only studied the behaviour of the two-scale numerical method in front of smooth initial 
data: it can be interesting to study the behaviour of the method on non-smooth initial data. 
Concerning the periodicity, we have built our numerical method in a periodic framework 
for mathematical reasons. However, it must be possible to extend this method to almost 
periodic framework in order to study problems like the noise generated by the blades of a 
turbine or a fan. If we go further, it can be interesting to improve the method in order to 
use it for studying the noise generated by multidimensional structures like cooling systems 
for computers or air conditioning, even if the links between time scales and space scales 
are not clear in this framework. 
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